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1 Introduction — astrophysical jets 

High velocity highly collimated beams of plasma - the jets - are known as a general 
phenomenon among astrophysical sonrces of different energy and spatial scales. At 
the npper end of the scale are the extragalactic jets emanating from active galactic 
nnclei. Little brothers of extragalactic jets has been detected recently as Galactic 
snperlnminal sonrces - the so-called microqnasars |l[]. These jets originate in stellar 
mass black hole accretion disks in high-energy binary systems. Protostellar jets are 
non-relativistic with velocities of several lOOkm/s [|^. 

The idea of a common astrophysical jet formation scenario is motivated by the ob¬ 
servational fact that jet formation is always connected to (i) the existence of (strong) 
magnetic fields and (ii) the presence of an accretion disk. The magnetohydrodynamic 
(MHD) model of jet formation (cf. ^) allows to describe the varions scales of 
astrophysical jets by a scaling of the jet magnetization. The theory tells us that 
highly relativistic jets must be also highly magnetized and that rotating MHD 
flows are subject of a self-collimation process P, 0. In this paper, we ask whether the 
ultra-relativistic motion and strong collimation indicated for GRBs is achieved by a 
mechanism similar to the “classical” astrophysical jets. 

2 Magnetohydro dynamic jets 

Highly relativistic MHD jets originate in the accretion disk surrounding a central 
black hole (Fig. 1). The disk provides the Poynting flux and drives the electric current 
system. The jet is initiated as a slow wind from the inner disk by a process which 
is not yet completely understood, in particular its time-dependent character. Most 
probably, some disk instability is responsible for the ejection of knots perpendicular 
to the disk surface. The disk wind is hrst launched magneto-centrifugally 0 and 
further accelerated and collimated into a narrow beam by Lorentz forces. The parallel 
component Fl,\\ ~ j± x is accelerating (or decelerating), while the perpendicular 
force Pl,± ~ j\\ x B is the collimating component (Fig. 1). Theoretical modeling. 
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Figure 1: {Left) The model of MHD jets launched from an accretion disk surrounding 
the central mass. {Right) The accelerating and collimating Lorentz force components 
are projected with respect to the magnetic flux surfaces T(i?, Z). 


therefore, requires to solve the governing MHD equations. In this contribution, we will 
concentrate on solutions to the stationary, axisymmetric, ideal MHD equations in the 
relativistic limit. The advantage of the stationary model is the possibility to obtain 
global solutions for the jet structure on spatial scales and resolution which cannot 
(yet) be reached by time-dependent simulations [§, Q. Long-term MHD simulations 
can, however, demonstrate the self-collimation of MHD jets (Fig.2; ||^). 


2.1 Relativistic stationary MHD 

The axisymmetric stationary ideal MHD equations can be reduced into two governing 
equations - the Grad-Shafranov equation (GSE) describing the force-balance across 
the magnetic flux surfaces^ T(r, 9) and the wind equation (WE) considering the forces 
along the held. Both equations are interrelated, however, in the relativistic limit of 
a high magnetization the held structure is almost force-free and can be investigated 
independently of the how dynamics. Here, we will discuss solutions of the force-free 
GSE in Kerr metric appropriate for jets from rotating black holes. 
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^ The magnetic flux is 9) = (I/Stt) f Bp ■ dA. In general relativity, the integration is over a 
loop of the Killing vector ffi = uZVf. In Minkowski space the flux is integrated over a circular area 
around the symmetry axis. 
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Figure 2: Non-relativistic MHD simulations of jet formation. Density (colors) and 
poloidal magnetic field lines after 0, 50, 100, 200 and 400 inner disk rotations. 


(see 1^). In Eq. (1), u) is the angular velocity of a zero angular momentum observer 
(ZAMO), a the lapse of ZAMO proper time to the global time, a) the cylindrical 
radius, and d)L(4/) the position of the two light “cylinders” (he.). Df( 4/) can be 
considered as rotational velocity of magnetic field lines. /(4/) is the total poloidal 
electric current. The general relativistic wind equation for the poloidal velocity u-p = 
Fup/c isQ [|, 0,0 

^ k^k2 - 2k2M^ - k^M^ 

{ko-M^y ’ ^ ^ 

with the Alfven Mach number M. Specihe total energy and angular momentum 
density E(\h), T(\h), are conserved quantities along T. As a key parameter for MHD 
jets, the magnetization (7(4/) measures the Poynting flux in terms of mass flux. 



a = 


4Mc3’ 


( 3 ) 


with the magnetic flux function $ = BpR^ and the mass flux M. The launching 
of a relativistic jet requires at least one of three conditions - rapid rotation, strong 
field, or low mass load (Eq. 3). In general, relativistic jets are highly magnetized 
|T3| , 1^, 0. Fixing the held distribution, one hnds a relation between the asymptotic 
Lorentz factor Fqo and a. We have (Michel-scaling) Foo = ~ (D| <F^/Mjet)^^^ for 

a spherical how with negligible gas pressure (“cold”) [0. In the case of collimation, 


the power law index is diherent. If the magnetic hux surfaces open up faster than 
spherically radial, the asymptotic flow is dominated by kinetic energy [O . 


^ with the following abbreviations, fco = 533^1 + + goo, k 2 = 1 — flplL/E), and = 

— (533 + ^goflL/E) + gooiL/Efl) / (^Qg - googso) 
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Figure 3: Solution of the relativistic MHD wind equation. Poloidal velocity u-p 
and particle density n over Goldreich-Julian density ncj versus radius. At the fast 
magnetosonic point, R = 33i?L5 both solution branches intersect. 


3 MHD jets in GRBs — structure and dynamics 

3.1 Parameter estimates 


We now estimate MHD jet parameters in the framework of GRBs assuming that these 
jets originate close to a stellar mass black hole and reach Lorentz factors up to 1000. 


Disk magnetic field: Differential rotation amplihes the toroidal held up to equi- 
partition, < Rgq- A disk dynamo may amplify poloidal helds up to Bp < B^. In 
advection dominated disks, Rgq — lO^G (M/5MQ)“^/^(Ma/ME)^^^(-R/3rg)“^/'‘, with 
central mass M, accretion rate Ma in terms of Eddington Me, and the horizon at 
15|. GRB models considering lO^^G helds are not consistent with this estimate. 

Light cylinder: The l.c. of the jet magnetosphere is located at Rp ~ c/VLp ~ 
4 X 10®cm (M/SMq) if the jet originates just outside the marginally stable orbit. 


Magnetization: The upper limit magnetization can be derived assuming the jet 
origin at the marginally stable orbit (highest held strength, most rapid disk rotation). 
Jet mass how rates ~ 10“^^MQyr“^ are indicated by GRB afterglows constraining 
the baryonic “jet” mass to IO^^Mq. This gives a jet magnetization of about 10^. 


3.2 Acceleration — asymptotic Lorentz factor 

Figure 2 shows a solution of the cold relativistic MHD WE in the Minkowski limit^ 
{a = 1000, Z{R) ~ The how accelerates rapidly, becomes super Alfvenic 

(R > Rp), super fast magnetosonic (R > 33Rl) and reaches asymptotically up ~ 300. 
The density remains above the Goldreich-Julian density, conhrming the underlying 
MHD assumption a posteriori. For cold jets, a is a free parameter. Thus, we may 
investigate the jet dynamics for diherent magnetization (Fig. 3). The two chosen hux 
distributions deviate from a spherical how and indicate the modihed Michel-scaling. 

^ For hot wind solutions considering the influence of Kerr metric on the jet acceleration, see 11 
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Figure 4: (Left) Modified Michel-scaling, Mp,oo(o') for BpR^ ~ R and ~ 

R~^'^(bottom). (Right) Global magnetic field structure of the collimation jet ^(R, Z). 


3.3 Collimation 


Figured shows a solution of the force-free GSE extending from the inner he. close to 
the hole, across the outer he. to the asymptotic, cylindrically collimated jet (a = 0.8, 
4/max = lO^^Gcm^, /max = lO^^A, i?L = lO^cm). The shape of the collimating jet is 
determined by the he. regularity condition. The held structure close to the hole and 
disk indicates a hollow jet with the bulk of the mass how in the outer jet layers E 


3.4 MHD jets in GRBs — open questions 

The numerical solutions presented here show that the classical MHD jet formation 
scenario can indeed be extended to GRB parameter in order to achieve both ultra- 
relativistic velocities and strong collimation. However, serious questions remain un¬ 
answered, among them: (i) Due to the Michel-scaling, it is difficult to account for a 
wide range in Lorentz factors, 10^ < F < 10®, which seems to be needed to generate 
the temporal variability of GRBs. (ii) Also, the bimodal distribution of GRB duration 
does not ht into the single MHD jet model. It is not clear, how such questions can 
be answered within a MHD jet model or whether pure electromagnetic or plasma 
processes are actually dominant. 
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Discussion 

Blackman: Can you comment on the claims of Ustyugova et ah that the boundary 
shape strongly influences the presence or absence of collimation? 

Fendt: We did not encounter this problem, maybe due to different boundary condi¬ 
tions in ZEUS code. Certainly, Ustyugova et ah’s point should be carefully considered. 

Hujeirat: Do you predict the existence of a super-Keplerian layer above the disk, 
especially as you said you do agree with the Blandford & Payne 1982 scenario? 

Fendt: The sub-Alfvenic disk wind co-rotates with the magnetic held. The toroidal 
velocity increases linearly ^ R along the how. Beyond the Alfven point, ~ R~^. 
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